Probing the central engine of long gamma-ray bursts and hypernovae with 

gravitational waves and neutrinos 

Yudai SuweE 

Department of Physics, School of Science, The University of Tokyo, 
Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan and 
Max- Planck- Institut fiir Astrophysik, Karl-Schwarzshild-Str. 1, D-85741 Garching, Germany 

Kohta MurascQ 

Yukawa Institute for Theoretical Physics, Kyoto University, 
Oiwake-cho, Kitashirakawa, Sakyo-ku, Kyoto, 606-8502, Japan 
CK ■ (Dated: December 3, 2009) 

O 

(<~*) ' There are the two common candidates as the viable energy source for the central engine of long 

, gamma-ray bursts (GRBs) and hypernovae (HNe), neutrino annihilation and magnetic fields. We 

investigate gravitational wave (GW) emission accompanied by these two mechanisms. Especially, we 
focus on GW signals produced by neutrinos from a hyper-accreting disk around a massive black hole. 
We show that neutrino-induced GWs are detectable for ~1 Mpc events by LISA and ~ 100 Mpc 
by DECIGO/BBO, if the central engine is powered by neutrinos. The GW signals depend on the 

• viewing angle and they are anti-correlated with neutrino ones. But, simultaneous neutrino detections 
are also expected, and helpful for diagnosing the explosion mechanism when later electromagnetic 
observations enable us to identify the source. GW and neutrino observations are potentially useful 

■ for probing choked jets that do not produce prompt emission, as well as successful jets. Even in 
non-detection cases, observations of GWs and neutrinos could lead to profitable implications for the 
central engine of GRBs and HNe. 

^ ' 

Qh' PACS numbers: 04.30.-w, 14.60.Lm, 98.70.Rz 

o 

liHj; I. INTRODUCTION 

, Gamma-ray bursts (GRBs) are the most luminous explosions in the universe. The production of GRBs is believed 

• to require that only the small amount of matter is accelerated up to ultrarelativistic speeds and collimated as a jet 
i 01- The duration of GRBs ranges from ~ 10"'^ s to ~ 10'^ s, vifith a roughly bimodal distribution of long GRBs 

OO ' of T > 2 s and short GRBs of T_< 2 s. The geometrically corrected gamma-ray energy of long GRBs is typically 

CO ■ ^ -^^^0-^2 |-,j,gg (ggg a^jgQ which is much smaller than the isotropic gamma-ray energy E]^° ^ ergs. 

\^ , The requirements for the energy budget and time scales suggest that long GRBs involve the formation of a black hole 

' (BH) (e.g., [^) or magnetars via a catastrophic stellar collapse event. The discovery of supcrnovae (SNe) associated 

0^ . with GRBs brought us the more direct evidence that GRBs result from a small fraction of massive stars that undergo 

■ a catastrophic energy release event towards the end of their evolution 0, [l] . Interestingly, some of the SNe associated 
J> \ with GRBs (e.g., GRB 980425, GRB 060218) showed evidence for broad hues indicating high-velocity ejecta with 

' inferred explosion energy of i?kin ^ 10^^ ergs (e.g., [^), and those energetic SNe are often called hypernovae (HNe). 
rS \ The central engine of both GRBs and HNe has to supply such an enormous explosion energy (for reviews, e.g., [lo|). 
^ . For the class of long GRBs, the commonly discussed candidates are massive stars whose core collapses to a black 
■ " " ' hole, either directly or after a brief accretion episode, possibly in the course of merging with a companion. This 
collapsar scenario is one of the most widely believed scenarios to explain the huge release of energy in GRBs and 
HNe [1, [m, ■ In this scenario, the collapsed iron core of a massive star forms a temporary disk around a few solar 
mass BH and accretes at a high rate (~ 0.1 — 10 Mq s~^), which is believed to produce a powerful jet leading to a 
GRB. The duration of the burst in this scenario can be related to the fall-back time of matter to form a disk or the 
accretion time of the disk. Possibly, HNe that are brighten by nickel could also be explained by a disk wind, which is 
subrelativistic with a speed comparable to the escape velocity of the inner disk p^ . 

Provided a disk and BH form, the greatest uncertainty in the collapsar scenario is the mechanism for converting 
the disk binding energy or BH rotation energy into directed relativistic outflows. So far, two general mechanisms 
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have been proposed: neutrino annihilation and magnctohydrodynamical (MHD) mechanisms in various kinds. In 
the former case, neutrino pairs are generated in the hot disk and impact one another with the largest angles along 
the rotational axis and deposit some fraction of the accretion energy [l^, [3, [3, [13, [3 [3 ■ It should be noted 
that the efficiency of energy deposition is no greater than ~ 1% of the total neutrino emission [2^, so that the 
required neutrino energy is E^, ~ 2.0^^"^^ ergs in order to achieve the jet energy of Ej ^ 10^^ ergs (which should be 
larger than the gamma-ray radiation energy). One important property of this dense, hot accretion flow is that the 
density is so high that cooling of the flow is dominated by neutrinos. This accretion flow is called neutrino-dominated 
accretion flow (NDAF). Although neutrinos can escape the flow more easily than photons, they can also be trapped 
in the flow and fail to escape before being swallowed by the black hole when the accretion rate becomes very high 
0, H, [3l- The other possible mechanism is a MHD-driven explosion, where two options exist for the energy source: 
the binding energy of the accretion disk [U, [l^] and the rotational energy of the BH through the Blandford & Znajek 
mechanism [2^, H^, H^- In the former case, the magnetic fields are amplified up to ~ 10^^~^^ G by rotation 
via some disk instabilities, leading to a strong outflow. In the latter case, the BH and accretion disk are connected 
with magnetic field lines and the rotational energy of the BH is extracted via these magnetic field lines. In some 
versions of the MHD scenario, the extracted energy propagates as the Poynting flux and is probably converted into a 
matter-dominated form at later stages, so that the efhciency of energy deposition is expected to be higher than that 
of the neutrino-annihilation mechanism. 

Although the various possibilities including the fast rotating magnetar scenario [13, [H, [2^, [sO, [Slj as well as the 
coUapsar scenario have been proposed, there has been no direct evidence for the central engine of GRBs. In this paper, 
we consider gravitational waves (GWs) as a probe of the jet-producing mechanism of GRBs (see [3l,[3l] for GWs from 
jet itself). The difficulty in probing the central engine comes from the fact that most observed properties of GRBs 
come from electromagnetic signals produced in regions far away from the central engine. On the other hand, GWs and 
neutrinos can reach us without losing information of physical conditions near the central engine. Especially, we focus 
on the GW emission generated by anisotropic neutrino emission, and estimate its amplitude which were investigated 
in the context of massive-star collapse [13, [H, [H, [13, [H, [11] (see also, [13, [HI)- The amplitude of this GW emission 
is governed by the radiated neutrino energy that has to be large enough for the neutrino-annihilation mechanism to 
work as the central engine of GRB jets. As a result, we can show that neutrino-induced GWs are detectable for ~1 
Mpc events by LISA and ~ 100 Mpc by DECIGO, if the central engine is powered by neutrinos. The coincident 
neutrino signals can also be detected by using the Supcr-Kamiokandc detector or recently proposed Mton neutrino 
detectors such as Hyper-Kamiokande. Possible successful detections combined with later electromagnetic observations 
can give us important clues to the central engine of long GRBs. 

This paper is organized as follows. In §11, we describe the formulation used in this work. The necessary expressions 
of the GW amplitude is presented. In §111, we give spectra of GWs from NDAF and apply the derived formula to 
a burst or multiple bursts representing GRBs. We also discuss detectability of GWs and whether we can identify a 
GW event as the formation of NDAF making a GRB and/or HN. In §IV, the gravitational wave background is also 
discussed. 



II. GRAVITATIONAL WAVES FROM NDAF 

Here we consider the gravitational wave emission from aniso trop ic neutrino radiation, which is proposed by Epstein 
[4^ . Following the formulation developed by previous works |34l . [ssl [sst , we first derive useful formulae of the GW 
amplitude for axisymmetric emission of neutrinos. 

First, we introduce the two coordinate frames shown in figure [TJ the source coordinate frame {x',y',z') and the 
observer coordinate frame {x, y, z). In these coordinate frames, the z'-axis is chosen as the symmetry axis of source, 
while z-axis is set to a line-of-sight direction. Furthermore, the origin of these two coordinate frames are set to 
the center of a GRB engine. Thus, the viewing angle denoted by ^ is given by the angle between the z-axis and 
the z'-axis. For simplicity, we assume that the z-axis lies in (a;', z')-plane. In this case, the two polarization states 
of GWs satisfying the transverse-traceless conditions become = h'^'J = —hyj and = K^y = K^J in the 
observer coordinate. The geometrical setup shown in Figure [T] yields the following relationship between the two polar 
coordinate frames {0,(f)) and {d',<j>'): 

sin 61' cos = sin 6* cos cos ^ + cos 6* sin ^, 
sin 0' sin 0' = sin sin 0, 

cos 6'' = — sin 6* cos sin ^ 4- cos 6' cos ^. (1) 



FIG. 1: Source coordinate frame {x',y',z') and observer coordinate frame {x,y,z). The observer resides at a distant point 
on the 2-axis. The viewing angle is denoted by 5, which is the angle between z and z' axis. The z'-axis coincides with the 
symmetry axis of the source. 



With these coordinates, one can obtain the GW amphtude as [35| 

hl^{t,S,) = ^ f'^dt'j^ rff7'(l+cos0)cos(2</.)^^^i^i^, (2) 

where G is the gravitational constant, c is the speed of light, D is the distance between the observer and the source, 
and dLi,/dQ' represents the direction-dependent neutrino luminosity per unit of solid angle in the direction of ft'. 
Note that the counter part of the amplitude, /i^"^, is obtained by replacing cos(20) by sin(20) in Eq. which 
vanishes when dLi,/dfl' is axially symmetric. 

In the above equation, and (j) are required to be expressed in terms of the angles 6' and (/)' with respect to the 
source coordinate variables, and the viewing angle, ^. Using Eq. ([T]), one can obtain the amplitude as following. 



2^ r^-D/c . dL^{n\ t') 



h+it,o^^l dt' / d^i'^{e',4>',i) \ (3) 



where ^'(6'', (/i',^) denotes the angle dependent factor, 

T ,/ >N /i ni ^ ■ ni ,1 ■ (sinS' cosit!)' cos£ — cos6'' sin£)^ — sin^ 0' sin^ (t!)' 

^-(0', </.',£) = (l + cos0'cos£ + sm0'cos0'smO7-— j-, ) . ^ , . 2 ■ 2 j^r W 

(sm&' cos (p' cos t, — cos t'' sm 4 ) + sm &' sm 0' 



This is a generalization of the Eqs. (26) and (27) of Miiller & Janka [35|, who derived the formula for the cases 
£ = and tt/2, respectively. By integration with respect to the azimuthal angle (</>') assuming that dL^{S},' , t')/dVl' is 
axisymmetric, we obtain 

h+{t,o = ^j dt' sine'de' ^{e',o^{9',t'), (5) 

where 

-27r[l + cos6''(2 + cos£)]tan2 (- ) (for 6 > ^) 
-27r [1 + cos e' {-2 + cos ^)] cot^ M (for g < ^ 
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When ^ = 7r/2, the above equations coincide with Eq. (8) of Kotake et al. [38[. With these equations, we estimate 
the amphtude of gravitational waves. For practical evaluations, we need to know dL^{^l,t)/dfl which should be 
determined by the nature of the central engine or produced jets. In the following subsections, we regard the shape of 
NDAF as a geometrically thin disk or an oblate spheroid, by which we obtain the angular dependence of dLi,{n, t)/dVl 
explicitly. 



A. A geometrically thin disk model 

We first consider a geometrically infinitely thin disk as a simple but useful case in order to model NDAF. In 
this case, the angular dependence of the neutrino emission is written as dL^/dfl oc \cos9\ if we assume that the 
emission of neutrinos is isotropic at the disk surface. We can write the neutrino luminosity per solid angle as 
dL^/dO' = [L^/2'k) I cos 0' I, where is the total neutrino luminosity. 

At first, we consider the case of an observer located in the equatorial plane for simplicity, which means ^ = 7r/2. In 
this case Eq.© is written as 



OC rt-D/c pTT IT 



-r- \ dt'LJt') / d6''(-l + 2|cos6''|)sin6''|cos6''| 
2G /1\ f'-^/" 



L,[t')dt'. (7) 



It is worthwhile to note that the anisotropy parameter a defined by Miiller and Janka [S^ (see Eq.(29) in their paper) 
is constant and equal to 1/3 in this case. Since the duration of neutrino emission is finite, the GW amplitude converges 
to a non- vanishing value [68|, which we denote h^Q. Witli cliRrRctcristic vrIucs, wg Ccld. writG hoo B^s 

where = L^{t')dt' is the total energy emitted by neutrinos. In this work, we employ E^y = 10^'' ergs for 

calculations of the GW amplitude. This is because the required energy for a GRB jet is Ej ~ 10^^ ergs and 
the efficiency of energy conversion is expected to be '--^ 1% for the neutrino-annihilation mechanism [l3l |. In addition, 
the energy released by matter at innermost stable circular orbit (ISCO) is about 40% of its rest-mass energy for an 
extremely rotating Kerr BH. This corresponds to ^ 10^''' ergs for accretion of 1 Mq. This value is also consistent 
with the observational evidence of GRB-HNe association because the explosion energy of HNe is 10 times larger 
(E'kin ^ 10^^ ergs) than that of ordinary core-collapse SNe, which emit E^ ^ 3x 10^^ ergs by neutrinos. // the explosion 
mechanism of HNe is similar to that of ordinary SNe, in which the explosion energy is provided by absorption and 
scattering of neutrinos, the total energy emitted by neutrinos must be larger than that of ordinary SNe, and we may 
expect El, ~ 10^** ergs. 

It is useful to see dependence of the GW amplitude on the viewing angle, ^. In the thin disk model, Eq.® can be 
integrated analytically, and we have 

2G l-h2cos£ o /£ 



f^-iO-^^^tan^^^) E.. (9) 

It is obvious that the GW amplitude is the largest when the observer is located in the equatorial plane, whereas 
the GW vanishes when the observer is located at the pole. This implies important implication, i.e., anti-correlation 
between GW signals and prompt gamma-ray photons. Prompt photons themselves are highly beamed (~ 1/r) since 
they are emitted from ultrarclativistic jets moving with F ~ 10^'^. Therefore, it may not be easy to detect prompt 
photons when we can expect the strong GW emission. Furthermore, there may be the possible population of failed 
GRBs whose jets do not generate prompt photons since, e.g., jets are choked in the stellar envelope. As discussed in 
the next section, GW signals would still be useful as a probe of the central engine even in such cases. 



B. An oblate spheroid 



We now consider an oblate spheroid in order to mimic the thickness of NDAF. In this case, the neutrino luminosity 
per solid angle is dLy/dVt cx \/ a? + {1 — a?) cos^ 6*, where a (< 1) is the ratio between the length of the major and the 
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minor axes. The gravitational wave amplitude for an observer in the equatorial plane = 7r/2) can be written as 



2G Vl -a2(l +a + Aa^) - 3a^{l + a)F{a) 




(10) 
(11) 



Fig[2]shows the GW amplitude from spheroid as a function of a with D ^ 10 Mpc and = 10^ ergs. For spherical 
emission (a = 1), no GW is are emitted by neutrinos. The case of the disk height being similar to the disk radius 
corresponds to a 0.5 that leads a suppression of GW by a factor of 3 compared to the case a = 0. 

In Fig. [3l we show hao{S,,a) normalized to its value for ^ = 7r/2 and a = 0, which is calculated numerically. The 
amplitude for ^ = tt/A and a = 0.5 is ^ 0.18hoc{TT/2, 0). It should be noted that the suppression is not so large even if 
we include both the viewing angle and thickness of the disk, as long as we consider a < 0.5 and ^ > 7r/4, respectively. 
It should be noted that the amplitude of GW is anti-correlated with neutrino. The apparent luminosity of neutrino 
shows the largest value for the observer at pole. When we consider the geometrically thin disk that was discussed in 
previous subsection, neutrinos can not be observed for the observer at the equatorial plane but GW is observable. As 
for the oblate spheroid, neutrinos are observable from any direction. The detectability of neutrinos will be discussed 
in the next section. 



o 




FIG. 2; The GW amplitude from NDAF using an oblate spheroid model as a function of the ratio between the length of the 
major and the minor axes, a. We employ D = 10 Mpc and E,y = 10*'' ergs in this plot. For spherical emission (a = 1), no GW 
is are emitted by neutrinos. The case of the disk height being similar to the disk radius corresponds to a ~ 0.5 that leads a 
suppression of GW by a factor of 3 compared to the case a = 0. 



III. THE GRAVITATIONAL WAVE SPECTRUM AND APPLICATION TO BURSTS 

In this section, we first obtain useful formulae of the GW spectrum from NDAF using the expressions derived in 
the previous section. We now consider axisymmetric sources, so that the only non- vanishing component is As a 
result, the local energy flux of GWs can be written as 



dAdt 16nG 



where dA = D^dQ. is the surface element. This is the general expression for the GW flux. 

We have calculated the angular dependence of dL^{^,t) / dVl in the previous subsections. Then, we can obtain GW 
spectra for NDAF mimicked by an infinitely thin disk or an oblate spheroid for given L^(t). As an example, let us 
obtain the useful expression for the case of a geometrically thin disk. (Numerical calculations are required for the 
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case of an oblate spheroid.) Let us insert Eq. ([7]) taking into account the angular dependence by Eq. ©. Then, we 
get 



dAdt 367rc5i:>2v y2 
The total energy flowing through dA between t = — oo and t — +00 is therefore 



l + 2cose)'tanM^ L,(t)^ (13) 



'^^^'^ ^ q + 2cosC)2tan4f^ M dtL,{t)\ (14) 



dA iGnc^D^^ ^' \2 

Integrating over a sphere surrounding the source, we find the total energy emitted by GW as 

1 ^ POO 

Egw - / dtL,{tf, (15) 

where /3 = ^ — 201n2 0.47039. This is consistent with Eq. (31) in Miiller & Janka [1^ since the anisotropy 
parameter a = 1/3 in this case as already noted in §11 [t^I ■ In order to get a GW spectrum, we write L^, in Eq. pS]) 
in terms of the inverse Fourier transform as 

/oo 
Z.(/)e-2"-^*d/. (16) 
-00 

After several algebraic steps, we easily obtain the GW energy spectrum as 

L.{f)\ (17) 



dEo^if) _2 G 



df g'^cs 

for a geometrically thin disk model for NDAF. 

Once we have a GW spectrum dE^^ /df (obtained numerically or analytically), we can discuss detectability of the 
GW. For this purpose, we shall use the characteristic GW strain expressed as 



for a given frequency / [ijl- Since we have this characteristic strain hdj) for either model of NDAF (a thin disk or an 
oblate spheroid), we can also compute signal-to- noise ratios (SNRs) obtained from matched filtering in terrestrial/space 
gravitational wave experiments. The optimal (with respect to the relative orientation of a source and detector) SNR 
is given by 
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where hn{f) = [^fSh^f)]^^"^ is the noise ampHtude with Sh{f) being the spectral density of the strain noise in the 
detector at frequency /. The factor five refers to averaging over all orientations of the source [1^. This determination 
of the SNR is optimistic, because a complete set of wave templates required for the matched filtering analysis is not 
available for burst events. However, SNRs for the other techniques of analysis (e.g., bandpass filtering and noise 
monitoring in Ref. [1^) can be connected with the SNR obtained from matched filtering within some factors. Thus, 
Eq. (fT9|) may give a little bit optimistic, but not severely overestimated values. 

In the following subsections, we give explicit expressions of L^{t) in order to apply our formulae to bursts. First, we 
show GW spectra regarding a GRB as one single burst caused by NDAF, where neutrino emission with ^ constant 
is assumed to occur during the duration of the GRB. However, the realistic time dependence of Ly{t) should depend 
on the nature of the central engine, which is unknown. In fact, the observed complex variability may be one of the 
clues to the behavior of the central engine. For example, the internal shock model, which is one of the most widely 
believed scenarios for explaining prompt emission, requires that jets are highly inhomogeneous (i.e., the dispersion of 
the bulk Lorentz factor should be large). The origin of many subshells moving with different speeds may be or may 
not be related to the activity of the central engine. In the former case, we may expect that neutrino emission occurs 
intermittentlyfTlj. This motivates us to consider the case of multiple bursts. 



A. A single burst 

In this subsection, we consider a single burst event. We assume the following time evolution of the neutrino 
luminosity 

i.(0- Y^e(t)e(r-t), (20) 

where T is the duration of the central engine of a GRB and Q is the Heaviside step function. This expression 
corresponds to the case where the neutrino luminosity from NDAF only depends weakly on time during the burst. In 
this case, the characteristic strain evaluated by Eq. (fT5|) can be written as 

/»c(/) = ^|sin(.r/)|, (21) 

where hao is the converged value of the GW amplitude obtained by Eq. ([8|). The spectrum converges to a constant 
value of v^^oo/"' for / <C l/T. This behavior typically arises from a sudden change of the neutrino luminosity and 
is called zero- frequency limit of GW memory [44| . 

Figure |4] shows typical examples of the characteristic strain he- In this figure, the red solid line represents the GW 
from GRB, which erupts at the center of Galaxy. The green dashed and blue dotted lines show the GW from a source 
at a distance of 10 Mpc but different durations, T = 10 and 200 s, respectively. The gray dot-dashed line shows the 
GW from a source with the same parameters as the green dashed line but Ei, ~ 10^^ ergs. Sensitivity lines for future 
planned interferometers, that is, advanced-LIGO, LISA, DECIGO/BBO, and ultimatc-DECIGO, are also shown in 
this figure [t^] . One can roughly read off the SNR from figure [H It should be noted that the GW amplitude is larger 
for the model with = 10^^ ergs and T = 10 s than one with E^, = 10^"^ ergs and T = 200 s in the range of / > 0.1 
Hz. Therefore, both the duration and the total energy emitted by neutrinos are important for the discussion of the 
detectability. 

More detailed information of the SNR can be found in figure [51 which gives the SNR for LISA. In this figure, the 
colored region corresponds to SNR> 1. The dotted Hue is for SNR^IO, which is required to detect the GW emission 
by the burst memory as discussed in Ref. [4^. As the duration becomes longer, the SNR of GWs from GRBs gets 
smaller. This is because the cutoff frequency, beyond which GW spectra depend on frequency as f~^, appears at 
/ ^ T^^ so that longer lasting events lead to smaller amplitudes at the high frequency range, whereas amplitudes 
in the low frequency range are independent of the duration. From this figure, one can see that LISA can detect 
GW signals up to distances of Z? ^ 1 Mpc. The same analysis shows that GRBs with distances of ~ 100 Mpc are 
detectable with DECIGO/BBO as long as the duration T < 10 s (see figure E]). 
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FIG. 4: The characteristic amplitude of GWs from a GRB. Thick lines represents characteristic amplitudes of GW for typical 
models. The red solid and green dashed lines are for T =10 s but D = 10 kpc and 10 Mpc, respectively. The blue dotted line 
is for D = 10 Mpc and T=200 s. Here, we assume the total energy to be E^, = 10^* ergs for red solid, green dashed, and blue 
dotted lines and = 10^^ ergs for the gray dot-dashed line. 




FIG. 5: SNRs of GW from NDAF for LISA. The white region in the right bottom part of the figure shows SNR< 1. The 
dashed line corresponds to SNR=10. 
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FIG. 6: Same as Figure El but for DECIGO/BBO. 
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B. Multiple bursts 



In this subsection, we consider the case of multiple bursts in order to take into account the intermittent time 
variability of the central engine. We employ the following time evolution for the neutrino luminosity, 

M<)^|;|i;e(<-i,T)e(-ii- + «-*). (22, 

where N is the number of subbursts and 6t is the duration of one subburst. We assume St to be constant for any 
subburst. The total duration is T + ST ~ T, since we consider the case of ST <C T [tS]. Then, one finds the 
characteristic strain as 

hcif) 



TT^NStf 



siniirStf) sin(7rr/) 



sin(7rT//iV) 



(23) 



Figure [7] shows the GW spectrum of Eqs. ([2T|) and (|23|) . The red sohd line represents the single burst case, which is 
the same one as the green dashed line of figure [D The green dashed line in this figure shows the case of multiple bursts 
with St = 0.005 s and N = 200. Both of the spectra are shown for T = 10 s. We can see that those two spectra are 
different in the high-frequency range because multiple bursts are caused by many bursts with the shorter time scale. 
On the other hand, these spectra coincide with each other in the low-frequency range because the long-term behavior 
is independent from the detail of the burst. The difference typically appears only above / ~ N/T ^ 20 Hz. Noting 
that the most sensitive frequency of LISA is about 1 mHz, we see that the detectability by LISA (and DECIGO/BBO) 
depends mainly E^, and D, so that the results of low-frequency GWs described in the previous subsection are unaltered 
even if a time variation is considered. By contrast, the detectability by advanccd-LIGO is affected by fluctuations 
with short timescales. The SNR for advanced-LIGO calculated by Eq. HH) is SNR - 9(£'^/10'5'^erg)(i:)/lMpc)-i for 
St = 0.005 s and N = 200, which is about 100 times larger than that for the single burst case. As a result, the GW 
might be detectable by advanccd-LIGO depending on the time properties of short-term variations during the burst. 
This implies that it may be possible not only to detect GWs by the burst memory, but also to extract important 
information about the temporal behavior of the central engine through GWs. 




FIG. 7; GW spectra of a single burst (red solid line) and multiple bursts (green dashed line) from NDAF. 



C. Can we identify NDAF as a gravitational wave source? 

Clearly, it is very important to identify sources emitting GWs in order to probe the central engine of GRBs and HNe. 
However, when we expect strong GW emission, prompt (electromagnetic) emission may be off-axis emission which is 
significantly diminished compared to on-axis emission. When the jet angle 9j is small enough compared to the viewing 
angle ^, the peak flux scales as ieFs)^^^. oc [1 — (V/c) cos(^ — 6j)]~'^^^, where V is the velocity of the jet [i^. Then, the 
peak flux of the off-axis emission is very roughly estimated as (e-FDmax ~ 10~^ ergs cm~^ -^7max 51 5(^/10Mpc)~ 



10 



for ^ ~ 7r/4 and Oj ^ 0.1, where i^l^^ax 51 5 ^^'■^ luminosity around peak energy (Throughout this subsection we 
define Qa ~ Q/IO" for a quantity Q incgs unit). Although the details of the spectra and light curves depend on the 
jet structure, viewing angle and so on [46L l47l| . they may be detected by the BAT on Swift or the GBM on Fermi, if a 
bright GRB occurs within 10 Mpc. However, we r nay also keep in mind that many of the discovered nearby GRBs are 
dimmer than classical high- luminosity GRBs [H, [l^]. For example, if the luminosity of such a low-luminosity GRB 
(for on-axis observers) is dimmer than that of a high-luminosity GRB by more than three orders of magnitude, and 
if the jet is well coUimated [t^ , it may be difficult to detect prompt emission even for very nearby GRBs. 

If we can detect both of the GW and prompt emission, we may infer the distance to the source by follow-up 
observations [75| . and then we can obtain some important information on the energy budget of the central engine 
such as E,y and the efficiency of energy deposition from . Although we may expect detectable prompt X or gamma 
rays from such nearby GRBs whose GW signals are detectable, we can think of a different possibility, that GWs and 
neutrinos are expected, but no prompt photons. In the coUapsar scenario, the jets have to penetrate the progenitor 
star to make the prompt emission. But, if they are chocked inside the progenitor, we do not expect any prompt 
emission (which is called a failed GRBs) [s^l- The possible population of such failed GRBs may not be neglected, 
as their local rate might be much higher than the rate of classical high-luminosity GRBs. Even in such cases, we 
may still expect detectable GW emission as long as the jets are launched by the neutrino-annihilation mechanism. 
In addition, coincident neutrino emission from NDAF should be expected, if GWs are caused by the neutrinos from 
NDAF. They may potentially help us to pin down the position of the bursts, and give us a clue whether GWs are 
caused by anisotropic neutrino emission (provided we know the distance). In figure [51 we plot detectabilities of 15-35 
MeV neutrinos by the Super-Kamiokande (SK) and Mton detectors as well as the SNR of GW signals expected by 
LISA. Pi is the cumulative Poisson probability to detect at least one neutrino event, X]^i(A*"6~^/n!), where fi is 
an expected event number, and P2 is the cumulative Poisson probability to detect more than one neutrino event, 
^^2(m"6~^/'^')i ill the 15-35 McV range. Thick lines refer to Mton neutrino detectors, while thin lines refer to 
22.5 kton neutrino detectors such as SK (see [5l| for the case of core-collapse SNe). For simplicity, we have used 
a Fermi-Dirac distribution for neutrino spectra with a temperature of 3 MeV, assuming that the neutrino radiation 
energy is shared between electron and anti-electron neutrinos. As we can see, coincident neutrinos can also be 
detected. Especially, if we have Mton detectors, we should detect coincident neutrinos if the GW emission is caused 
by neutrinos, and possible detections of scattering events (rather than capturing events shown in the figure) become 
useful for positioning (see [s^, for the discussion of simultaneous observation of neutrino and GW in the case of 
core-collapse SNe). 
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FIG. 8: The detection probability of neutrinos and SNR of GWs from NDAF. We assume E„ = 10^'' ergs. See text for detail. 

Later (follow-up) observations of photons seem to be crucially important to determine the distance to the source, 
as has been already mentioned. First, observations in the optical and radio bands may lead to discovering a HN (e.g., 
see [13]) associated with the formation of NDAF making jets. If GRB jets are successful and make prompt emission, 
possible detections of the prompt photons would be helpful for this purpose. In addition, it may also be possible to 
detect afterglow emission (which can be expected not only for successful GRBs, but also failed GRBs [111). Whether 
the outflow is well coUimated or not, it should be decelerated by colliding with the interstellar medium. When prompt 
emission is not observed, afterglows arc called orphan afterglows. In this case, the afterglow is also dim in the early 
phase before the jet break, during which it may not be so easy to detect its emission. However, the expected flux 
after the time when T{tp)£^ — 1 docs not depend on whether afterglows are orphan or not [H, [5^. Following the 



11 



external shock scenario @], we can roughly estimate the peak flux after the peak time tp ^ 8.1 x 10^ s [5 + 21n(10) + 
21n(e/0,,-i)](C/0,,-i)'</'i(i;ir53)'^''^"'/'> as (eF,)_, ~ 6 x 10"^ ergs cm'^ s"! £^,-1 E^^°,,{D/Mpcyh^' for 
KT53 ~ 10^^ ergs, - 0.1 and 71 - 1 cm , where is the fraction of nonthermal electron energy, E^j is the kinetic 
energy of the ejecta, and n is the density of interstellar medium [5^, respectively. Although actual detectabilities 
depend on the detail of the afterglow spectra, this afterglow emission could be detected by groundbased optical 
telescopes and X-ray satellites such as Swift, future MAXI and EXIST. Since the spectra and light curves depend 
on 9j and ^, their detection could potentially give us some information about relevant parameters such as 6j and 
^, leading to a crude estimate of the asphericity parameter a under our NDAF scenario. Those later observations 
(i.e., detecting an associated HN component, prompt and afterglow emissions), combined with GW and neutrino 
observations, may help us to reveal the central engine of GRBs and HNe. 



IV. THE GRAVITATIONAL WAVE BACKGROUND 



We are now in a position to discuss the contribution of GWs from NDAF of GRBs to the background GW radiation. 
According to Phinney [131 , the sum of the energy densities radiated by a large number of independent GRBs at each 
redshift is given by the density parameter flcwif) = Pc^i^PGW /dlog f) as 



^Gw(/) 



1 



PcC^ Jo 



dz 



l + z 



dt 
dz 



df ' 



(24) 



where pc = 3H'^/{8nG) is the cosmic critical density, i?GRB (z) is the GRB rate per comoving volume, and fz = (l+z)/. 
The cosmological model enters with \dt/dz\ = [(1 + z)H{z)]~^ and, for a flat geometry, H{z) = HqI^Ia + ^Im{^ + z)^]- 
We adopt the cosmological parameters f2j\f = 0.3, JIa ~ 0.7, and Hq — 100 ho km s~^ Mpc~^ with hg = 0.72. For 
the GRB rate, we employ the following GRB rate history (GRBS model in Refs. [H, [5§|) in units of yr~^ Gpc~'^, 



-Rqrb — Rq- 



46e 



3.4z 



,3.8z 



45 



(25) 



where F(z, flm, ^a) = \/^a + ^^m(l + z)^/{l + z)^!"^ . Note that the dependence on the overall GRB rate history is 
so small that it is enough to use the above history for our purpose of demonstration. 

In figure [9l the calculated r^cw is plotted together with the sensitivity curves of future detectors. We assume 
that all GRBs have T = 10 sec for simplicity. The red solid and green dashed lines are obtained for a GRB rate of 
i?o = 18 which is expected for high-luminosity GRBs, while the blue dotted line is with one for i?o = 500 which may 
be expected for low-luminosity GRBs. The red solid line is obtained with the single-burst spectrum (Eq. (1H|)), and 
the green dashed line is for the multiple-burst spectrum (Eq. ([23]) ) so that these lines differ above 1 Hz. The blue 
dotted line is for the multiple burst spectrum. The enhancement of the GW background due to inclusion of time 
variation does not change the detectability by space-based interferometers since the most sensitive frequency is about 
0.1 Hz. 

Although we have used the inferred rates for successful GRBs, we may expect a larger contribution to the GW 
background if there is large population of failed GRBs. At present it is difficult to predict the rate of such choked 
jets, but we might expect that the rate of failed GRB is typically smaller than the SN Ibc rate. Then, the possible 
contribution to the GW background is at most JIqw^o ^ 10^^^ at ^ 0.01 — 0.1 Hz (see [5^ [6^ for a discussion of 
the neutrino background from GRBs and HNe). 



V. SUMMARY AND DISCUSSIONS 



In this paper, we have calculated the amplitudes and spectra of GWs from gamma-ray bursts, which are expected 
in the neutrino-annihilation mechanism above an NDAF disk in the coUapsar scenario. We show that GWs can be 
detected by LISA and DECIGO/BBO if the source distance is Z? < a few Mpc. The expected rates of high- and 
low-luminosity GRBs occurring at _D < 10 Mpc, which correspond to the rates of bursts that can be observed by 
^ 5 Mton neutrino detectors [54I, are ^ 2 x 10~^ yr~^ and ~ lO^'' — 10""^ yr~^, respectively. Hence, detections of 
GW signals are possible only when GRBs fortunately occur in the nearby universe. Nevertheless, the possible GW 
signals are very important in the sense that they can probe the mechanism of the GRB central engine. Obtained 
GW spectra are fiat at low frequencies, which is the characteristic feature of the GW memory effect. In order to 
detect GWs with LISA or DECIGO/BBO, sufficiently large values of the anisotropy parameter a and large values 
of the released burst energy are typically required. Since the MHD mechanism may typically lead to a smaller burst 
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FIG. 9: The energy density parameter of gravitational wave background (GWB). See text for detail. 



energy compared to the neutrino-annihilation mechanism (where E^, ^ iQ^^-b-i ^^^^ needed), detections of flat GW 
spectra at low frequencies from extragalactic GRBs would allow us to infer that the central engine is governed by 
the neutrino-annihilation mechanism above an NDAF disk. We show that coincident neutrinos should be detected if 
we have Mton detectors as long as GWs arise from the neutrino emission, which should also give us important clues 
about the central engine (e.g., estimate of Ei, allows us to make a consistency check for the burst energy triggering 
the GW emission). 

In order to identify the source emitting GWs, later electromagnetic observations will be crucial. We may see prompt 
emission, afterglow emission and supernova emission. Although the former is significantly diminished when the viewing 
angle is very large, we may be able to detect it for very nearby bursts. Afterglow emission is also detectable even if we 
have orphan afterglows. Successful observations of GWs, neutrinos and photons will allow us to diagnose the central 
engine of GRBs and to obtain some clues on the jet production mechanism. In addition, we can expect that HNe 
are associated with GRBs (but it may not be necessarily the case). Detections of HNe in the optical and/or radio 
band will give us profound information about the explosion mechanism of HNe. Strong GW and neutrino emission 
can be expected since the neutrino-annihilation mechanism requires that the radiated neutrino energy is very large 
{E,^ ~ 2^q53-54 gj-gs). On the other hand, it would be more difficult to expect detectable GW signals from the MHD 
mechanism. Hence, if we have very nearby GRBs, non-detections of both GWs and neutrinos might suggest the MHD 
mechanism rather than neutrino-annihilation which seems to have some caveats (e.g., the very strong dependence on 
the disk temperature and inefficiency in production of fireballs). Whether signals are detected or not, simultaneous 
observations by GWs and neutrinos will be important to probe the central engine of GRBs and HNe, even though 
the burst rate is not so large. 

In this work, we have used a thin disk or an oblate spheroid model to mimic NDAF. Such a simplified treatment 
is sufficient for our purpose, and the characteristic feature in the low-frequency range will not be altered even if 
we employ elaborate numerical simulations. But, although sufficiently realistic calculations have not been done so 
far, further studies and developments may allow to evaluate both the neutrino flux and the GW amplitude more 
quantitatively. Also, the effect of neutrino oscillation, which is neglected in this work, should be included for more 
quantitative predictions. In addition, we have employed the fiducial energy emitted by neutrinos as E^, = 10^^ ergs, 
which is also unknown. The assumptions employed for this value are that the energy of a GRB jet is Ej ~ 10^^ ergs 
and the efficiency of energy conversion from neutrinos to a jet is 1% for the neutrino-annihilation mechanism. Since 
GRBs have a variety of energetics, E^ might be significantly larger or smaller than this value. Detections of GWs 
and neutrinos are useful since the signal-to- noise ratios of GW show a linear increase with Ei^. In addition, we might 
have information about the efficiency of energy conversion, which is theoretically important, through obtaining the 
jet energy from electromagnetic observations. 

Although we have mainly considered production of GRB jets with a duration of ~ 10— 100 s, required for explaining 
prompt emission, there are additional possibilities of the GW emission. Recent observations by Swift have shown that 
many GRBs have flares in the afterglow phase, which suggest that the central engine is active for a longer time than 
the duration of the prompt emission [6l| . Flares are also likely to be energetic, and the radiation energy of some flares 
is even comparable to that of the prompt emission. Possibly, they may be produced by slower and less collimated jets, 
and their geometrically corrected energy may also be comparable to that of the prompt emission. Then, we could 
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also expect GW emission even in the early afterglow phase. 
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